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We theoretically consider infrared-driven hyperbolic metamaterials able to spatially filtering terahertz radia- 
tion. The metamaterial is a slab made of alternating semiconductor and dielectric layers whose homogenized 
uniaxial response, at terahertz frequencies, shows principal permittivities of different signs. The gap provided 
by metamaterial hyperbolic dispersion allows the slab to stop spatial frequencies within a bandwidth tunable 
by changing the infrared radiation intensity. We numerically prove the device functionality by resorting to full 
wave simulation coupled to the dynamics of charge carries photoexcited by infrared radiation in semiconductor 
layers. 



Manipulating terahertz (THz) radiation is generally a 
difficult task since the most of standard materials simply 
do not respond to such frequencies. However, the advent 
of metamaterials has allowed to partially reduce this dif- 
ficulty since their electromagnetic properties can be arti- 
ficially manipulated 1 through a suitable design of the un- 
derlying constituent unit cells. At the same time a num- 
ber of setups have been proposed for steering the THz 
radiation 2 and reconfigurable electrically 3 or optically 4 ^ 
driven metamaterials have been exploited for conceiving 
active THz devices^. 

The most of the proposed active THz devices are tun- 
able frequency-domain filters since it is relatively sim- 
ple to control the metamaterial dispersion properties 
through external stimuli. In this Letter we theoretically 
propose a way for achieving active and spatial filtering 
of the THz radiation by means of a suitable hyperbolic 
metamaterial whose THz response can be tuned by an 
auxiliary infrared field. Hyperbolic or indefinite media 7 
are uniaxially anisotropic metamaterials having princi- 
pal permittivities of different signs, a remarkable feature 
leading extraordinary plane waves to be ruled by a hy- 
perbolic dispersion relation. Hyperbolicity is the main 
physical ingredient leading to unusual optical effects as 
negative refraction 8 and hyperlensing 9 and supporting a 
number of proposed devices as beam splitters 10 , spatial 11 
and angular filters 12 and optical switches 13 . The tun- 
able hyperbolic metamaterial we consider in the present 
Letter, together with the fields geometry, is sketched in 
Fig.l. The metamaterial slab of thickness L is obtained 
by stacking along the x-axis alternating layers of an in- 
trinsic semiconductor 14 (sc) and a negative dielectric (nd) 
of thicknesses d sc and d n d, respectively and it is illumi- 
nated by an infrared (IR) plane wave linearly polarized 
along the ?/-axis and normally impinging onto the slab 
interface at z = 0. The THz field (TH) is a transverse 




FIG. 1. (Color online) Layered metamaterial setup together 
with THz (TH) and infrared (IR) plane waves geometry. 



magnetic (TM or p-polarized) monochromatic plane im- 
pinging with incidence angle onto the interface. 

The infrared field within the semiconductor layers pho- 
toexcites electrons to the conduction band which dynam- 
ically recombine so that the resulting electron density N 
is described by the rate equation^ 5 - 



(1) 



where H is the Planck constant divided by 27r, eo is the 
absolute vacuum permittivity, € sc (ujir) is the semicon- 
ductor permittivity at the infrared frequency ujir, E^ ir ^ 
is the infrared field within the semiconductor layers, tr 
is the carriers' nonradiative recombination time and B 
is the coefficient of radiative recombination. The semi- 
conductor dielectric behavior at the infrared frequency is 
described by the linearized permittivity model 15 
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where uir is the infrared semiconductor refractive index 
background, No is the transparency value of the carrier 
density, Sno the refractive index change at transparency, 
c is the speed of light in vacuum and A is the differen- 
tial absorption coefficient. At steady state (dN/dt = 0), 
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FIG. 2. (Color online) Terahertz homogenized dielectric per- 
mittivities of Eqs.© as functions of the normalized local op- 
tical intensity of the infrared field. 




Eqs.flTJ) and © yield 
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where W = \E^ IR ) \ 2 /\E sat \ 2 and \E sat \ = 
[rReofiiRcA/ (2Hujir)] 1 , so that electron density shows 
a saturable behavior (i.e. TV ~ N if > \E aat \), 

as a consequence of the saturable absorption model of 
Eq.(j2]). Maxwell equations for the field E^ IIV) together 
with Eqs.(j2j) and (j3j) describe the infrared nonlinear 
behavior within semiconductor layers. The resulting 
(and self-consistently evaluated) electron density N has 
a strong impact on the semiconductor response at the 
THz frequency ujth, since the permittivity is 

6 SC (U TH ) = 6^(UJ TH ) + - ' (4) 

eo uth{± - iujtht) 

where e s °J (uoth) is the semiconductor permittivity in the 
absence of the infrared radiation, — e and m* are elec- 
tron charge and reduced mass and r is the relaxation 
time. Note that the term proportional to N in Eq.(j4]) is 
the standard Drude permittivity contribution due to con- 
duction band electrons and it is the physical ingredient 
allowing the infrared field to tune the overall slab THz 
response. If the layers thicknesses d sc and d nd are much 
smaller than the THz wavelength (of the order of tens of 
microns), the overall metamaterial slab of Fig.l shows a 
homogeneous THz uniaxial response with permittivities 
(the optical axis lying along the staking x-axis) 16 
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FIG. 3. (Color online) a) Reflectivity R and b) transmissivity 
T at the wavelength Xth — 23.08 \ivn as functions of the 
THz incidence angle for a homogenized slab of thickness 
L — 9.23/xra and permittivities of Fig.l at three different local 
infrared intensities I^ IR \ c) Terahertz hyperbolic dispersion 
curves for the three cases of panels a) and b). The hyperbola 
vertex is at k x — ky — (27t/Ath)a 



j(iR) — €oC \e( ir } | 2 /2 (normalized with the saturation 
intensity I sat = eoc\E sat \ 2 /2 = 3.7 kW/ cm 2 ), with the 
simplified assumption that such intensity can be regarded 
as uniform within the bulk of the slab and equal to inci- 
dent infrared plane wave intensity (see below). The con- 
sidered wavelengths are Xir = 2ttc/ujir = 0.879 jim and 
Xth = 2ttc/ujth — 23.08 \im and we have chosen gallium 
arsenide (GaAs) as semiconductor (for which tr = 10 ns, 
B = 1.3 • 10- 10 cm 3 s-\ n IR = 3.6, 5n = -0.07, N = 



9.91 -10 17 cm" 3 , A= 1.71 -10- 15 cm 2 , €^(cj T h) = 12.9, 
r = 3.29 • 10" 13 8 and m* = 0.067m , m being the elec- 
tron mass) with filling fraction f sc = 0.1 and a negative 
dielectric with permittivity c n d(wTH) — —0.61 + 0.003i 
(this value coinciding with that of calcium fluoride CaF2 
with imaginary part reduced by a tenth for discussion 
purposes). It is evident from Fig. 2 that the structure 
has been tailored to show indefinite permittivity tensor 
in the absence of infrared illumination, i.e. Ke(e x ) < 
and Re(e^) > for j( Ji? ) = 0. Note that the homoge- 
nized structure keeps its indefinite character for increas- 
ing /( /jR ) and the permittivities globally show saturation 
at high intensities as a consequence of the carrier density 
saturation predicted by Eq. ([3|) . Remarkably, the infrared 
intensity /( /jR ) has a strong impact on Re(e^). In Fig. 3 we 



have plotted the reflectivity R = E^ 7 ^ / 



E 



(TH) 



and 



fndtnd^TH) (5) the transmissivity T 



where f sc = d sc /(d sc + d nd ) and f nd = d nd /(d sc + d nd ) 
are the layers filling fractions whereas e n d{oJTH) is the the 
dielectric permittivity of the negative dielectric layers for 
which Re [e nd (ujTH)} < 0. 

It is well known that, since the permittivities of Eqs.(j5j) 
involve averages of different kinds, it is possible to tailor 
the structure in such a way that e x and e z have differ- 
ent signs, the resulting hyperbolic medium begin tun- 
able, in the present analysis, through the auxiliary in- 
frared field. In Fig. 2 we have plotted e x and e z evaluated 
from Eqs.(|5]) as functions of the infrared optical intensity 



(TH) 
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(TH) 



(see Fig.l for 



the definition of the field amplitudes) as functions of the 
THz incidence angle for a homogenized slab of thickness 
L = 9.23 fim with permittivities coinciding with those 
of Fig.2 at the intensities iW = 0, iW = I sat and 
j( IR ) — 10J sa t. It is worth noting that all the three trans- 
missivity curves have a bell-shaped profile with maximum 
located at an angle dependent on the infrared intensity 
(analogously the reflectivity have a complementary be- 
havior) so that the slab can be regarded as an active spa- 
tial filter allowing (forbidding) transmission (reflection) 
of those THz plane waves with incidence angle close to a 
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FIG. 4. (Color online) THz transmissivity R and reflectivity 
T of the metamaterial of Fig.l, evaluated trough full waves 
simulations, for three different values of the optical intensity 



T (IR) 



of the incident infrared plane wave. 



central angle in turn tunable through the infrared optical 
intensity. 

The physical mechanism supporting such an active 
spatial filtering functionality can easily be grasped by 
noting that the chosen TM incident THz plane waves 
couples to the slab extraordinary waves whose disper- 
sion relation is k%/\e z \ — k%/\e x \ = (2tt/ Xth) 2 , where the 
signs of the permittivities (from Fig. 2) have been explic- 
itly reported for clarity purposes (neglecting permittiv- 
ities' imaginary parts). Such dispersion relation is, in 
the k x k z plane, a hyperbola with vertices along the k x - 
axis located at k x — ±ky where ky = (2tt/\th) \/\^z\ 
(see panel (c) of Fig. 3 where the hyperbola correspond- 
ing to the considered three infrared intensities are plot- 
ted). Momentum match at the interface z = implies 
that the extraordinary plane waves have transverse wave 
vector k x = (2tt/\th) sin# so that in the two situations 
|sin0| < \Z\ez~\ or |sin0| > \f\e^\ (i.e. \k x \ < ky or 
I fee | > ky : respectively) the externally impinging THz 
plane wave excites, within the slab, evanescent or propa- 
gating waves, respectively. As a consequence, the trans- 
mission is low for sin < ^/\e z \ (where evanescent waves 
provides a residual "tunnelling" radiation), it reaches a 
maximum at greater angles (due to the slab propagating 
waves) and it eventually vanishes at 6 = 90° for geomet- 
rical reasons. 

The discussion of the THz spatial filtering function- 
ality has hitherto been based, for clarity purposes, on 
three simplified assumption i.e 1) low loss regime (small 
imaginary part of the negative dielectric permittivity), 
2) THz electromagnetic homogenization of the layered 
medium (validity of Eqs.flSJ)) and 3) uniformity of the 
infrared optical intensity within the slab bulk. While 
the first two assumptions can easily be supported, the 
third one is more serious since the layers periodicity is 
generally comparable with the infrared wavelength and, 
in addition, reflection of infrared radiation by the slab 
together with its nonlinear behavior within semiconduc- 
tor layers have to be considered. In order to account 
for of all these physical ingredients and to show that the 
structure actually behaves as an active THz spatial fil- 



ter, we have resorted to full wave simulations where lin- 
ear and nonlinear Maxwell equations for the THz and 
the infrared field are coupled to the electron dynamics 
described by Eqs.(|3]) and We have chosen the layers 
thicknesses d sc = 24 nm and d n d = 220 nm and the per- 
mittivities e n d((jOTH) = —0.61 + 0.03i, € n d(0JlR) = 2.04 
which are the CaF2 permittivities at the considered THz 
and infrared frequencies, respectively, whereas all the re- 
maining involved parameters are those used above. The 
results of the simulations are reported in Fig. 4 where 
we plotted the THz transmissivity and reflectivity of the 
structure for three different values of the optical inten- 
sity l\ IIV> = eoc|E^ /jR ^ | 2 /2 of the incident infrared plane 
wave (see Fig.l). Note that the overall spatial filtering 
functionality is evidently exhibited by the considered re- 
alistic structure, even though the values of the trans- 
missivity are smaller as compared to those of Fig. 3 as a 
consequence of the CaF2 absorption. 

In conclusion we have proposed a metamaterial struc- 
ture which, driven by an auxiliary infrared field, is able 
to spatially filtering THz radiation. We believe that the 
proposed combination of hyperbolic metamaterial and 
semiconductor concepts can suggest different ideas for 
conceiving novel and efficient THz active devices. 
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